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In this paper a real option approach for the valuation of real assets is 
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presented. Two continuous time models used for valuation are described: 
O geometric Brownian motion model and interest rate model. The valuation 

for electricity spread option under Vasicek interest model is placed and the 
formulas for parameter estimators are calculated. The theoretical part is 
confronted with real data from electricity market. 

\Q ; PACS numbers: 02.30.Cj, 02.50.Ey, 02.70.-c 

oo 

Q . 1. Introduction 

^ ' The liberalization of electricity market caused that modeling on this mar- 

ket became very important skill. It helps us to minimize loss and hedge our 
position. It is a very interesting fact that spread option could be used for 
valuation some real assets as power plants or transmission lines. But before 
that we need to know the spread option price formula. Very popular model 
used for option valuation is the geometric Brownian motion model but it 
is not very efficient. In the 2005 the idea of modeling domestic electric- 
ity market using interest rate model was introduced by Hinz, Grafenstein, 
^ ' Verschue and Wilhelm pQ. They valuated European call option written on 

power forward contract under Heath Jarrow Morton model and, in this way, 
created very interesting class of models. 

The aim of my work is the valuation and calibration of electricity spread 
option under interest rate model applied for electricity market. I start with 
assumption of Vasicek model and using martingale methodology |2| valu- 
ate spread option. Using the maximum likelihood function methodology I 
estimate model parameters. I compare constructed model with geometric 
Brownian motion model by applying both models to real option valuation. 
I make simulations to show the difference between to discussed models. 



(1) 



2 



EBS printed on February 2, 2008 



My paper is organized in the following way. At the beginning (Section 2) 
I describe what does it mean that we bought a spread option, next (Section 
3) I introduce the reader into real options world. In Section 4 I describe val- 
uation methodology for spread option under interest rate model and present 
also option price formula for geometric Brownian motion model. The cal- 
ibration methods for both models are described in Section 5. At the end, 
in Section 6, all theoretical deliberation are confronted with real data and 
some simulation results are presented. 

2. Electricity Spread Options 

In this section two interesting cross commodity derivatives on electricity 
market are described. The first one is the spark spread option, which is based 
on fact that some power plants convert gas into electricity. The underlying 
instrument is the difference between the gas and electricity prices (the spark 
spread). The basic parameter connected with this kind of instrument is the 
heat rate, the ratio which describes the amount of gas required to generate 
1 MWh of electricity. The definition of such an instrument has a form E] : 
An European spark spread call option written on fuel F, at fixed swap ratio 
K, gives its holder the right, but not the obligation to pay K times the unit 
price of fuel F at the options maturity T and receive the price of one unit 
of electricity. 

It is easy to imagine such kind of option which better fits the Polish 
electricity market. We should assume that the underlying instrument is the 
difference between the carbon and electricity prices. But in this time there 
is no possibility of valuation of such an option because we don't have the 
representative carbon price. Generally, if we assume that Pe and P F are 
respectively future price of lMWh of electricity and the future price of the 
unit of fuel and K is the swap ratio than we could describe the payoff of the 
European electricity-fuel spread call option as 

C f (Pe,Pf,T) = max[P E (T) - KP F (T),0] 

and the payoff of the European electricity-fuel spread put option has form 

P F (P E ,P F ,T) = max[KP F (T) - P E (T),0] 

The second derivative is the locational spread option. It is based on 
fact that transmission of power from one location to another is very popular 
transaction. It is normal, for transmission system, that the power is moved 
from the place of lower price to the place of higher price and this is why the 
transaction is profitable. The whole transaction depends on the difference 
between the electricity prices and also on delivery costs and for hedging we 
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could use options. This kind of instrument could be denned in following way 
[3]: An European call option on the locational spread between the location 
one and location two, with maturity T, gives its holder the right but not the 
obligation to pay the price of one unit of electricity at location one at time T 
and receive the price of K units of electricity at location two. Assume that 
Pi and i"2 are the electricity prices at the first location and second location 
respectively. The payoff of the European locational spread call option is 
given by 

C L (Pi ,P 2 ,T) = max[Pi (T) - KP 2 (T) , 0] 

The put option is defined similar and the payoff of the European electricity- 
fuel spread put option has form 

P L (P 1 ,P 2 ,T) = ma X [KP 2 (T) - Pi(T),0] 



3. Real Options 

Suppose that for describing the commodity we use three qualities [G, t, L): 
G - the nature of good, t - time when it is available, L - location where it 
is available. We could define a real option as technology to physically 
convert one or more input commodities (G, t, L) into an output commodity 
(G',t',L'). For example most of power plants are real option because they 
give us the right to convert fuel into electricity. The transmission line is also 
real option. It gives us the right to change the electricity in one location into 
electricity in second location. The works of Deng, Johnson and Sogomonian 
[HUS] contain two formulas defining how to valuate generation and transmis- 
sion assets. If we define that up(t) is a one unit of the time-t right to use 
generation asset we could say that it is the value of just maturing, time-t 
call option on the spread between electricity and fuel prices Cp(t) and the 
one unit value of capacity of power plant using some fuel F is given by 

V F = I u F (t)dt = [ C F (t)dt 
Jo Jo 

where T is the length of power plant life. 

Similar if we define that u AB (t) is a one unit the time-t right to convert 
one unit of electricity in location A into one unit of electricity in location B 
we could say that it is the value of just maturing,time-t call option on the 
spread between electricity prices in location A and B C£ B (t). The one unit 
value of such transmission asset is given by 

rj~i rj ~\ rj ~\ rj~t 

V L = [ ui B (t)dt+ [ ul A {t)dt= [ c£ B {t)dt+ [ CE A (t)dt 
Jo Jo Jo Jo 

where T is the length of transmission network life. 
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4. Valuation methods 

In this section I present the widely known geometric Brownian motion 
model and I valuate the call spread option for the new, interest rate model 
using martingale methodology. All calculations are described below. 

4-1- Geometric Brownian Motion Model 

Suppose that the future prices of commodity are described by following 
stochastic differential equations 

dP^t^) = fHP^t,^ + a 1 P 1 (t,T)dW tA , 

dP 2 (t, T) = fi 2 P 2 (t, T)dt + a 2 P 2 (t, T)dW t , 2 , 



where Wt t i = pWt )2 + \A ~~ P 2 W t2 and Wt t 2, W t2 are i.i.d. Brownian mo- 
tions. It is known fact p], that the price of the spread call option with swap 
ratio K and time to maturity T, written on futures contract with maturity 
U < T is given by 

C 1 (t) = e- r ^[P 1 (t,UMd + (t))-KP 2 (t,UMd^(t))}, 

where 

/„ Pi(t,U) - « 2 (T-t) 
d ± (t) = - in ^ 2 



aVT-t 
and 

2 2 2 

a = a i — 2a\pa 2 + a 2 . 

4-2. Interest Rate Model 

For domestic currency, for example MWh, we denote two processes: 
pi(t,T), p 2 (t,T) which are the future prices of one unit of commodity. The 
interest rate functions for such processes are respectively 

dr tt i = (ai - hr t) i)dt + o x dW t ± 

and 

dr t ,2 = (a-2 ~ b 2 r tj2 )dt + a 2 dW tt 2, 

where Wt,i = pWt t 2 + \A — P 2 W t 2 and Wt >2 , W t2 are i.i.d. Brownian mo- 
tions. We assume that there exist the savings security N t (for example in 
USD), with constant interest rate r, for which 

p {t T ) = p{t ' T) 
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is the USD price of future delivery of 1 unit of commodity. We know £Q 
that there exist a martingale measure P for which the discounted processes 
N a P2 n'^ are martingales. We have 

C 2 (0) = N E P (( Pl (T, U) - K P2 (T, C/))+iV- 1 |F )/e- rt iV 4 

We define the new discounting processes for 1=1,2 as 

dB t ,i = B tti r tti dt, 

where Bq^ = 1. If we change the measure from P to Pi 

cflP _ N T B ,i 
d¥l ~ N B T ,i 

we know that processes p~i(t, T) = P1 q^ , P~2(t, T) = P2 g t '^ are Pi-martingales. 
From interest rate theory we obtain 

dp~ 1 (t,T)=p~ 1 (t,T)n 1 (t,T)dWt,i, 

where 

ni(f,T) = -£(l-e-fcC r -*)). (1) 
If we change the measure again from Pi to P2 

dPi = B TA B 0t2 
dF 2 B 0tl B T , 2 

we know that p 2 (t,T) = and p \(t,T) = = p ~ l( ^ } f M are P 2 - 

martingales and similar to the earlier situation, we have 

dp 2 (t, T) = p 2 (t, T)n 2 (t, T)dW t , 2 , (2) 

where 

n 2 (t,T) = -^(l-e^ T -% (3) 

After simple calculations we also have 

dp\(t,T) =p 1 (t,T)n 1 (t,T)dW t>1 , (4) 

where W t)1 = P W t)2 + y/T^W^, and W' t>2 = W' t>2 + = t. We 

also assume that 

dN t = N t vdV t , (5) 
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where Vt = p-\Wt,i + \A — $Wt2-> an d W*,2, W" t ' 2 , W t " 2 are independent 
Wiener processes. For discounted processes following equation is true 

P[t T) _ Pi(t,T) _ pj(t,T) 

W)~ e -rt Nt ~ e _ rtjV W 

Having the necessary stochastic differential equations we could price the 
option. We change the measure from P 2 to Q in following way 

dP 2 = B T , 2P2 (0,U) 
dQ B , 2 p 2 (T,U)' 

Process X(t, T) = p^yj = is Q-martingale. From Ito lemma we 

know that 

dX(t,T) = X{t,T){n 1 {t,T)dW til -n 2 (t,T)dW t , 2 ), 

where W tl =pW t2 + JT^Wt 2 and W t2 = W' 2 + ^^C^'Mt, 

Now, for calculation of the option price, we could use the Black-Scholes 
formula 

p 2 ((WK° ff Pi(T,[7) n+|F ,_ p((] mp (( P2(T,U) + 

P 2 (0, U)Eq((X(T, U) - K) + \F a ) = P 2 (0, C/)(X(0, C/)$(d+) - tf$(d_)) = 
Pi(0, ?7)$(d + ) - KP 2 (0, U)*(d-), 

where 

d± = - * 2 



<j 2 (7,T) = y (n 2 (u,T) - 2n 2 (u,T)pni(u,T) + n|(u,T))di 



<7(0,?7) 

9/ m\ „ / m\ / m\ . 9/ m\\ , 
- 

It 

and $ is the normal cumulative distribution function. For every time point 
< t < T the option price with swap ratio K and time to maturity T, 
written on futures contract with maturity U < T is given by 

C 2 (t) = P ± (t, U)<Z>(d + (t)) - KP 2 (t, U)$(d-(t)), 

where 

, , , m KP 2 (t,U) * 2 
±( } = ^) ' 

This methodology could be used directly for locational spread options and 
also for fuel-electricity spread options if we assume that the swap ratio be- 
tween MWh and unit of fuel is one. 



EBS printed on February 2, 2008 



7 



5. Historical Calibration 

In this section we describe how to fit our models for real, historical data. 
At the beginning we assume that we are given historical prices of future 
contracts Pi(tk,Tj) and P^^^^Tj), k = 0, ...,n, j = 0,...,m, in discrete 
time points to < t\ < . . . < t n and To < T\ < . . . < T m , where t^+i —tk = dt 
and Tj + \ — Tj = AT . 

For geometric Brownian motion model the calibration methodology is 
not very complicated. We analyze the returns of future prices of instrument 
and n is its mean, a 2 is its variance and correlation parameter is simply the 
correlation between returns of two instruments. But for interest rate model 
the calibration is quite complicated, especially for multidimensional HJM 
model 0. Calibration for discussed Vasicek model is presented below. 
Let us consider following process 

n (f rp\ _ Pi(t,Tj) PjfrTj) 

m[hlj) Pi(t,T j+ x) m,T j+1 y 

From Ito lemma we know that 

dr)i(t,Tj) = r]i(t,Tj)[(n 2 (t,T j+1 ) - n i (t,T j+1 )n i (t,T j ))dt+ 

+(m{t,Tj) - ni(t,T j+1 ))dW t ,i] 

We could write that rii(t, T) = nj(T — t) because this function depends only 
from the difference between the maturity time T and the time point t. If we 
then consider the process 

' 1i{t,T) 

we know that Sj(T — t) is normally distributed with mean 

ai(T -t) = (n 2 (T + AT - t) - m{T + AT - t)m(T - t))dt 
and variance 

f3 2 (T -t) = (m(T - t) - m(T + AT - t)fdt. 
Knowing the form of functions rii{T — t) 0,® we see that 

(3 2 {T-t) = ^ e - b ^[l-e- biAT ]\ 2 dt (7) 
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and 



P?(T-t + dt) 

After discretisation and for assumption that dt = AT = 1, T — t = pAT 
and j = 1, . . . , m we could say that the estimator of j3 has the form 

.. m 

tfipAT) = - V7 S ? .(pAT) - Si j(pAT)si), 
j'=l 

where 

j m 

«i = — y~! s j,i(pAr) 

m ^— ' 

and we put 

PijTj-pAT,^) P i {T j -{p+\)AT,T j ) 



Sij(pAT) = r,iTj ; ' A/V/ '' 1 ''"W/'+iiAT.^+A-n 



So using 0,(0) we have 



P l {T ] -pAT ) T j ) 
Pi (Tj —pAT,Tj+AT) 



1 A 2 (PAT) 



2AT A((p + 1)AT) 



A 2 (pAT)6, 

<7i ~ 



It is easy to notice that the correlation parameter between processes sx(T—t) 
and 52(7" — t) is p, so we have 

£™ iKifrAT) - ^(^-(pAT) - S - 2 ) 



E?=i(*ij(pAr) - JiJYEiLiC-jy (pAT) - s 2 ) 2 

At the end we should calculate also parameters connected with process N t . 
From equation (JJjJ we know that for i=l,2 

^ t ,T) = i^p- = e- rt m,T). 

Using Ito lemma and formulas @, Q, J21 we could calculate following 
dynamic 
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T) = T)[w 2 - unaCt, T)pi]dt + &(t, T)[(n 2 (i, T) - v Pl )dW t , 2 - 



We know that the process 



v(T-t) 



is normally distributed with mean (v 2 —vri2(t, T)p\)dt and variance (nf(i, T)- 
2vri2(t,T)pi + i; 2 )ctt so we have that for (it = 1 

v = n,2 (pATj + 2y 



and 



-2 - 

w - 2/ 

Pi 



where 

v ( v at) = e ~ rATp ^ zl^LIA - p ^ - (p + i)Ar,rj) 

y ' J KP ' Pi{Tj -(p+ 1) AT, Tj)AT 

and 

y = — Y]yj(pAT). 



6. Simulation and Conclusion 

For simulation I used data from New York Mercantile Exchange (NYMEX). 
I considered historical quotation of future natural gas (Henry Hub) and 
electricity (PJM) contracts since January, 2004 until March, 2006. The pa- 
rameters were calculated using calibration methods described before. All 
estimated parameters are presented in Table El I assumed that the constant 
interest rate is r = 0.05. For valuation of gas fired power plant I assumed 
that the life-time of the power plant is T = 15 years and Pe,o = 55.750 
USD, P Ffi = 6.3080 USD. 

In Figure 1. we see the value of power plant for the heat rate ranging from 
5 to 15 for both presented models. We could notice that there is difference in 
changes dynamic for analyzed models. The value of power plant for interest 
rate model is much more smaller than for GBM model and it tends to zero 
when the heat rate goes up. It is a very good feature, because in reality the 
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TABLE I 



Estimated parameters for GBM model and for interest rate model using 
historical data from New York Mercantile Exchange. 



Geometric Brownian Motion 



°9 
Pe 
Pg 
P 



1.0945 
1.2943 
4.4098 
4.8145 
0.8688 



Interest Rate Model 



0"e 


0.0678 


a 9 


0.0042 


be 


3.7515 


b 9 


1.8205 


P 


0.1892 


Pi 


0.7266 


V 


0.0668 
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Fig. 1. Top panel: Future prices of natural gas and electricity for contract maturing 
in March, 2006. Bottom panel: Simulated unit value of gas fired power plant, with 
life length 15 years, for both models. 



value of power plant for heat rate greater than rj 9 should be close to 
zero. Looking at work of Deng we could say, that the value of power plant 
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under GBM model is usually too high, so also in this aspect the interest rate 
model gives better results. 
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